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Abstract 

We study the Kozai dynamics affecting the orbital evolution of trans- 
neptunian objects being captured or not in MMR with Neptune. We provide 
energy level maps of the type (u, q) describing the possible orbital paths 
from Neptune up to semimajor axis of hundreds of AU. The dynamics for 
non resonant TNOs with perihelion distances, q, outside Neptune's orbit, a N , 
is quite different from the dynamics of TNOs with q < a N , already studied 
in previous works. While for the last case there are stable equilibrium points 
at u = 0°,90°, 180° and 270° in a wide range of orbital inclinations, for the 
former case it appears a family of stable equilibrium points only at a specific 
value of the orbital inclination, % ~ 62°, that we call critical inclination. We 
show this family of equilibrium points is generated by a mechanism analogue 
to which drives the dynamics of an artificial satellite perturbed by an oblate 
planet. The planetary system also generates an oscillation in the longitude 
of the perihelion of the TNOs with i ~ 46°, being Eris a paradigmatic case. 
We discuss how the resonant condition with Neptune modify the energy level 
curves and the location of equilibrium points. The asymmetric librations of 
resonances of the type 1:N generate a distortion in the energy level curves and 
in the resulting location of the equilibrium points in the phase space (co,q). 
We study the effect on the Kozai dynamics due to the diffusion process 
in a that occurs in the Scattered Disk. We show that a minimum orbital 
inclination is required to allow substantial variations in perihelion distances 
once the object is captured in MMR and that minimum inclination is greater 
for greater semimajor axis. 
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1. Introduction 



One of the open problems of the trans-neptunian region (TNR) is to 
find an explanation for the wide variety of orbits that the discovered trans- 
neptunian objects (TNOs) exhibit. In particular, those objects with per- 
ihelion outside Neptune's orbit and with high orbital inclinations and ec- 
centricities that cannot be explained by the diffusive process the Scattered 
Disk Objects (SD Os) experience, nor b y the past dyn amical history of the 



Gladman et al.l . 120021 ; iFernandez et al 



planetary system dTsiganis et al.l. 12005 ). It is known (IDuncan et al.l . [1995 



20041 ) that TNOs with perihelion 



q < 36 AU are inside a chaotic region that generates a diffusion in semi- 
major axis enhancing aphelion distances. But the diffusion itself cannot 
decouple the perihelion from Neptune's region, then it is necessary to in- 
voke another mechanism to explain the existence of high perihelion SDOs, 
also known as Detached Objects. The origin of these problem atic orbits have 
give n rise to theories involving passing stars (llda et al.l . l2000l ). scattered plan- 
ets (iGladman and Chanl. 120061) . tides from the star c luster where the Sun was 
formed ( Brasser et al. . 20061 ). a stellar companion ( Gomes et al. . 2006 ) and 
others. But, in order to avoid extra hypothesis it is necessary an abroad 
panorama of the dynamics generated by the planetary system itself in the 
TNR. Secular resonances, mean motion resonances (MMR) and the Kozai 
resonance (KR) are dynamical mechanisms operating in the TNR and that 
could explain some eccentric orbits decoupled from encou nters with Neptune 



([Fernandez et al.l . l2004t iGomes et al.l . 120051 ; iGomesl . 1201 ll ) . The present work 



contributes with new results following that line of thinking. 
Secular resonances do not affect the TNOs located beyond a 



42 (IKnezevic et al 



199ll ) and MMRs with Neptune do not generate themselves a substantial 
yariation in the o rbital elements, as for example the MMRs with Jupiter do 
( IGallardol . l2006bl ). The orbital variations observed when a particle is cap- 
tured in MMR with Neptune are due to a secular dynamics inside the MMR, 
like Kozai dynamics. The Kozai d ynamic s , or m ore properly Kozai-Lidov 
dynamics, has its roots in a study by iLidovl (119621 ) of the secular evolution of 
an artificial Earth's satellite perturbed by the Moon. Applying these ideas 
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to the asteroid belt. iKozail (119621 ) developed an analytical approximation for 
the mean secular perturbing function due to Jupiter assumed in circular or- 
bit, valid even for high inclination and high eccentricity asteroidal orbits. 
The conservation of the semimajor axis, as is usual in secular theories, but 
also the conservation of the parameter H = \/l — e 2 cosi, and the existence 
of an energy integral, K(ou,e), allowed the calculation of the energy level 
curves for this model showing large oscillations of e and i coupled with the 
argument of the perihelion, u, around equilibrium points that appear for 
H < 0.6 at oj — 90° and 270°. These maps of energy level curves have axial 
symmetry with respect to u = k90° being k = 0, 1,2,3,4 and with respect 
to i = 90°. The orbital regime under oscillations of u it is known since then 
as Kozai resonance (KR) and it provides an explanation to the existence of 
some asteroids in spite of their large orbital oscillations. 

The conservation of H gives us the ultimate limits in e and i between 
which the particle's orbit can evolve assuming a secular dynamical regime, 
that means no close encounters with the planets, whether resonant or not. 
But the actual limits in e and i for a specific orbit are imposed by the 
orbital variations allowed by the energy level curve to which the particle is 
confined. We show in Fig. [T]the known population of TNOs with q > 36 AU, 
that means outside the diffusion region, in a diagram H(i, e). These objects 
at present do not undergo encounters with Neptune and do not experience 
diffusion in semimajor axis, then their eccentric or very inclined orbits are due 
to the early history of the planetary system or to the secular dynamics acting 
at present. Each object could evolve along its corresponding if-constant 
curve allowing to increase or decrease its present perihelion distance, but 
the limits only can be defined when studying the energy level curve for each 
obj ect. 



Kozail (119851 ) modified his model to include also a mean motion resonant 



condition between the mean longitudes of Jupiter and the asteroid assuming 
a fixed value of the critical angle that characterizes the resonance, allowing 
to understand the KR inside a MMR (or MMR+KR). Energy level curves 
for resonant asteroids are very different than for the non resonant case, and 
also are different the limits for the orbital oscillations. Moreover, for some 
resonances, new stable equilibrium solutions appear for u different from the 
known 90° or 270°. These were called asymmetric equilibrium points some 
years after and the y appear as stability island s between collision trajectories 



with the planets (IGronchi and Milanil . Il999l ). Except for the case of the 



Jupiter's Trojans, the energy level diagrams exhibit, as in the non resonant 
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cas e, symmetry with r espect to u = k 90 c 



Bailey et al.l (119921 ) studied the orbital evolution of high inclination long 



period comets with a new semi-analytical averaging method that takes into 
account the perturbations of Jupiter assumed in circular orbit. The method 
is not based on any analytical series expansions but in calculating numerically 
the mean energy integral, doing it good for any set of orbital elements. They 
numerically computed energy level curves which long period comets should 
follow assuming no planetary close encounters or MMRs with Jupiter. That 
energy level curves showed stable equilibrium oscillations around u> = k90° 
for cometary perihelion inside the planetary region. The principal result was 
that some comets have extreme perihelion oscillations which is clearly related 
to the origin of some sungr a zing c omets. 



Thomas and Morbidellil ( 119961 ) followed an analogue method but includ- 



ing the four jovian planets in order to study the dynamics in the outer Solar 
System. The series of energy level curves they obtained give a very complete 
panorama of the perihelion behavior from the Nep tune region down t o the 
proximity of the Sun. They also obtained, as did iBailey et al.l ( 119921 ). the 
same equilibrium oscillations around u = k90° for the case of long period 
comets. Nevertheless, by their detailed study, it is clear now that the equi- 
librium islands at u = 90°, 270° with very small q are more stable than the 
ones at u = 0°, 180° which vanish in some circumstances due to the prox- 
imity with the collision curve with the planets, which is a curve in (u, e) 
that implies orbital i ntersection between the pa rticle and the planets for a 
gi ven value of (a.H). iMichel and Thomas! (119961 ) followed the same method 
of Thomas and Morbidelli ( 19961 ) but including the effects of all the planets 
allowing the study of the dynamics of asteroids with a < 2 AU in NEA type 
orbit, showing by first time the Ko z ai dyn amic s associated with the terres- 
trial planets. iKinoshita and Nakail ( 119991 ) and iKinoshita and Nakail (120071 ) 
finally provided analytical solutions for the original Kozai model of an aster- 
oid imd^r_aii_e2rtemal_pCTturber. 



Fernandez et al.l ( 12004 ) and I Gomes et al.l ( 120051 ) in a series of numerical 



simulations showed the relevance that the interaction between MMRs and KR 
have for increasing significantly the perihelion distances of SDOs, from q ~ 40 
AU to q ~ 70 AU. They found the KR generates large perihelion variations 
only when the TNO is cap tured in MMR wit h Neptune and preferably in 
resonances of the type 1:N. iGomes et al.l ( 120051 ) also provide figures showing 
ho w the energy level curves are modified when a MMR condition is imposed. 



Gallardol ( I2006af ) using a classical expansion up to order 10 of the dis- 
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turbing function due to Neptune analyzed how the KR is generated inside 
and outside a MMR with Neptune. It was found a new regime of the KR for 
orbits with perihelion outside Neptune in the far TNR: it appears a stable 
equilibrium point at u — 90° and 270° but for a very particular value of the 
inclination, i ~ 63°, called critical inclination, which re sembles the dynamica l 
behavior of Molniya like artificial satellites. Curiously, iKuchner et al.l (120021 ) 
had also detected notable eccentricity variations in fictitious particles with 
initial circular orbits with % = 61° at a = 43 and a = 45 AU but with u oscil- 
lating around 0° and 180°. Then, it is evident there is an interesting Kozai 
dynam ics related to orbits with % ~ 61° — 63° in t he TNR. IWan and Huang 
( 120071 ) using the same classical expansion used in iGallardol ( ]2006al ) due to 
Neptune studied the effects of the KR inside the MMRs 2:3 and 1:2 pro- 
viding energy level diagrams for that resonances, but they did n ot consider 
i nclina tions greater than 55°. We must take into account that iKrasinsky 
( 119721 ) proved the existence of the critical inclination in the framework of the 
restricted circular three-body problem for both situations: outer perturber 
and inner perturber, which are i ~ 39° and i ~ 63° respectively, for the case 
that the ratio of the semimajor axis of the inner body with respect to the 
outer body tends to zero. 

Just for completeness, the Kozai dynamics was also studied in extraso- 
lar systems as a natural extension of the original asteroidal problem, that 
means, the planet that generates the perturb a tion is outside the pe r turbe d 
planet, see for exam ple Innanen et al.l (1997 ). Wiegert and Holman ( 1997h. 



Hol man et al.l ( 119971 ). iMikkola and Tanikawal ( 119981 ) and lValtonen and Karttunen 



( 120061 . chapter 9), while the mutual perturbations of two-planet systems 

were only considered more recently in the context of the three body problem 

( Mi ght chenko et al.l . 120061 : iLibert and Tsiganisl . l2009t iMigaszewski and Gozdziewskil . 



2010l ). It is worth mention that a big deal of the analytical studies of the 



Kozai dynamics applied to satellites, small bodies, and exoplanets were de- 
veloped assuming an external perturber, being the case with a perturber 
inside the perturb ed orbit, as occurs in the TNR, only analyzed with de- 
tail more recently (IMigaszewski and Gozdziewskil . l2010t iFarago and Laskar , 



20101 : IMigaszewski and Gozdziewskil . l201ll ) . 

In the present work , follo wing the framework of the mo del given by 



Thomas and Morbidellil ( 119961 ) and lMichel and Thomas! (119961 ). we start ex- 



tending in Section [2] their study from Neptune region to large orbital a, q 
and i finding the possible limits of perihelion variations assuming a secular 
non resonant evolution and, especially, studying in detail the Kozai dynam- 
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ics around i ~ 62°. The results in form of energy level curves are presented 
in terms of (u, q) for an immediate comprehension of the possible perihelion 
variations. In Section [3] we present a secular analytical model that describes 
approximately the secular dynamics in the far solar system given an explana- 
tion for the oscillations around % ~ 62° and predicting new orbital behaviors. 
In Section H] we analyze how a MMR enhances the KR strongly affecting the 
energy level curves and, especially, the peculiar 1:N resonances with Neptune. 
In Section [5j analyzing a series of numerical integrations of fictitious particles 
with semimajor axes from 50 to 500 AU we discuss the applicability of the 
secular Kozai dynamics and the conditions for the capture in the MMR+KR 
and for the enhancement of the perihelion distances. 



2. Non resonant Kozai dynamics beyond Neptune 



We followed the method described in iThomas and Morbidellil ( 119961 ) but 
including all the planets from Mercury to Neptune assuming coplanar and 
circular orbits, as is usual. The inclusion of the terrestrial planets is not 
relevant for the study of the TNR but we have included them in order to 
study particular cases with q inside the inner planets if necessary. In this 
context we numerically calculated the double integral defined by the mean 
of the perturbing function over mean anomalies of the planets and the test 
particle and obtained a mean perturbing function which is independent of the 
particle's longitude of the ascending node, Q, and mean anomaly, M. This 
implies that a and H = \/l — e 2 cos i are constant as explained in the refer- 
ence above. We constructed energy level curves K{u j, q), where we preferred 



the pe rihelion distance instead of X = \/l — e 2 as in lThomas and Morbidelli 



( 119961 ) because it directly gives the path in q that a TNO can perform. 

In Figs. [2] and [3] we show a synthesis of energy level curves for H = 
0.1, 0.3, 0.5 and 0.7 and for a = 50, 100, 200 and 300 AU. They are symmet- 
ric with respect to u = k90°, and i = 90° so we limited the plot to the 
region 0° < u < 180° and i < 90°, or H > 0. As H ->• 1 (or e, i ->• 0) 
the energy level curves are progressively more straight horizontal lines indi- 
cating that the allowed variations in e, i become negligible. The structure 
for q < 30 AU showed by Fig s. [2] and [3] is well known since the work of 



Thomas and Morbidellil ( 119961 ). Just for illustration, we show in Fig. H] a 
typical map of energy level curves in the region interior to the giant planets. 
Although the origin of some centaurs could be explained with this dynamics, 
due to the energy level curves connecting q ~ with q < a N , we remark 
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that, analyzing all ours energy level diagrams that cover a wide range of q 
in the TNR, we conclude that, in general, there is no connection between 
q > 30 AU and q < 30 AU. That means, this secular, non resonant, dynam- 
ics does not allow stable perihelion variations that cross Neptune's region. 
There exist some trajectories that enter the region q < 30 AU but inevitably 
reach the collision curve with Neptune and in that case the method is not 
more valid, the energy is not more conserved and the motion will not follow 
the energy level curves. Nevertheless, we will see later in Section H] that the 
boundary q ~ 30 AU can be crossed in a smooth regular orbital evolution 
following the energy level curves in the case of TNOs cap tured in MMRs 
with Neptune as showed for example by iGomes et al.l ( 120051 ). 

We also note that in the region q > 30 AU in general there are not 
relevant variations in q except for objects with H < 0.5 where an impressive 
island of equilibrium appears at u = 90° and 270°, clearly shown in Fig. [3] 
in a narrow range of inclinations around 63°. By fine tuning the parameters 
(H, a) we reconstructed the transformations that these equilibrium islands 
experience from a = 36 to a = 500 AU which are showed in Fig. At the 
lower values of a for which the equilibrium points exist the corresponding 
inclination is ~ 61° and from a > 100 AU the inclination of the equilibrium 
point is ~ 63°. The perihelion variations that they generate can be of about 
10 AU or more, and always occurs for q outside Neptune. Fig. [5] shows that 
in fact the birth of the equilibrium points occurs for a ~ 36 AU and are 
located at i ~ 61°, u> = 0°, 180°. For larger a they appear new equilibrium 
islands at i ~ 61°, u = 90°, 270° which are the only ones that survive for 
a > 54 AU. Equilibrium points around u = 90° are present in previous works 
but for very different, lower inclinations and very low perihelion values, well 
inside the planetary region, contrary to the case showed in our Figs. E] and 
13 In Fig. [6] it is showed a fictitious particle evolving in the TNR obtained 
from a nu merical integration of t he outer Solar System covering 1 Gyr using 
EVORB ( iFernandez et al.l . 120021 ) and the corresponding energy level curves 
in a good match. Below the equilibrium points at i ~ 63° it is verified 
duj/dt > and above , doj/ dt < 0. These plo ts explain the previous results by 
Kuchner et al. fl2002h and lGallardol (l2006af ). It is evident that the oscillation 



of u, which defines the Kozai resonance (KR), is installed in all the TNR 
at i ~ 61° — 63°, constituting a family of equilibrium points in the Kozai 
dynamics. 

From our maps we can conclude there is a connection between perihelion 
near Neptune and very low perihelion, which is a well known result from 
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previous studies, and, more important, there is another connection between 
perihelion near Neptune and perihelion outside Neptune's orbit through a 
new equilibrium point located at i ~ 62°. We remark there is no connec- 
tion between both regions, that means, from the outer region to the inner 
one. Neptune's orbit is like a barrier between these regions. In fact, the 



chaotic evolution associated to orbits with q < 36 AU ( IDuncan et all 11995 



Gladman et all . l2002r Fernandez et all 12004 ) allows the connection between 



both regions, but not more under the hypothesis of a secular regular evo- 
lution. We will show in Section H] that this scheme is strongly modified for 
particles captured in MMR with Neptune, but first, in the next section we 
will find a dynamical justification for the equilibrium points at the critical 
inclination. 



3. An analytical secular non resonant approach 



We wi ll obta i n a sim ple secular model which is in essence the same pro- 
posed by iKozail ( 119621 ) but applied to the TNR and taking into account 
the perturbations due to N planets moving in circular and coplanar orbits. 
The secular evolution of a TNO with perihelion outside Neptune's orbit can 
be approximately modeled by the perturbative effects of a series of massive 
belts representing the N planets. For each planet j we define a belt in the 
plane (x,y) with radius equal to the semimajor axis of the planet, aj, and 
with mass equal to th e plane t's mass uij in units of solar masses. Following 
Brower and Clemencd (Il96ll chapter III) it is possible to show this model 
generates a perturbing function given by: 



C ix 



9E /j, 



35 z 4 



R = -r^ l - 3^ + —4 1 - lO^r + — -r + . . . 



4 r 



64 r 



3 r 4 



where C is the moment of inertia in the direction of z: 



(1) 



N 



C = J2 m i 



(2) 



and 



N 



m j a j 



(3) 



being \x = k 2 , the square of the Gaussian constant, and z and r are z com- 
ponent of the TNO and its heliocentric distance respectively. Actually, the 
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planets and also the Sun considered as a solid body contribute to C and E 
but the contribution due to the Sun is negligible compared with the contribu- 
tion of the planets, then we discard the Sun's contribution to the perturbing 
function. Taking into account the four major planets we obtain C = 0.1161 
M AU 2 and E = 52.0057 M AU 4 . The ratio between the term depending 
on E respect to the term depending on C is about ~ 300/r 2 which for a 
typical TNO is less than 0.1, but we cannot neglect it. Calculating the mean 
respect to mean anomaly of the perturbed TNO we get the mean perturbing 
function R m : 

R ™ = - 3 / 2 ,3/ 2 (^ + ^) + -.. (4) 
lba 6 1 — e z ) 



where 
and 

Ra = 



R 2 = C{1 + 3cos2i) (5) 
9E ((2 + 3e 2 ) (9 + 20 cos 2i + 35 cos 4i) + 40e 2 (5 + 7 cos 2i) cos 2u sin 2 i) 



512a 2 (1 - e 2 ) 2 

(6) 

being a,e,i the TNO's orbital elements. For a typical TNO, R4/R2 < 0.1 
and vanishes for grow i ng va lues of a. Using Laplace planetary equations 



([Murray and Dermottl . 119991 ) we obtain the secular evolution of the TNO: 



de 45ekE 2 

(5 + 7cos2i) sin i sin 2u + . . . (8) 



dt 512a n /2(i _e 2 )3 

di 45e 2 kE 
Jt ~ _ 1024a n /2(i-e 2 ) 4 

dVL 3Ck 



dt 4a 7 / 2 (l-e 2 ) 2 
doj 3Ck 



(5 + 7cos2i) sin2i sin2co' + . . . (9) 
cosz + ... (10) 



dt 16a 7 /2(i_e 2 ) 2 
and from Eqs. (fTUI) and (TTTT) : 

dm 3Ck 



(3 + 5cos2i) + ... (11) 



dt 16a 7 / 2 (l-e 2 ) 2 



(3-4cosi + 5cos2i) + ... (12) 
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where only the lower order terms are showed. Perturbations for 
are due to the lower order terms of the perturbing potential, that means, the 
terms depending on C, while the perturbations for i, e are due to higher order 
terms of the perturbing potential depending on E. Note that di/de — e(l — 
e 2 ) -1 / tan i which is the same relationship we can deduce from y/l — e 2 cos i = 
H, where H is constant. 

From these equations it follows that particles with orbits verifying cos(2i) ~ 
—3/5 (or % ~ 63° and 117°) will verify duj/dt ~ while, in general, de/dt ^ 
and di/dt ^ 0. In an energy level diagram (u,q) this evolution corresponds 
with a trajectory which has a vertical tangent (u> = constant) when that 
inclination is reached, indicating a returning point in the energy level curves. 
In particular, at u = k90° according to the equations above we also have 
de/dt = di/dt = 0, obtaining four equilibrium points, two stable and two 
unstable. In an alogy to the m otion of an artificial satellite there is a critical 
inclination, i c , (IBeutlerl . 120051 chapter 3) ( ICapderoul . 120051 . chapter 3) which 
produces dw/dt = 0. Above i c the argument of the perihelion has a retro- 
grade motion and below i r a prograde mot ion. This is in good agreement 
with the result o btained by lGallardol (l2006a| ) and with the value obtained by 
Krasinskyi ( 119721 ) using a different approach in the framework of the restricted 
circular three-body problem. 

Substituting e 2 by its equivalent 1 — (H/ cosi) 2 in Eq. (j3J) we obtain R m 
depending on two variables (u,i) and two parameters (a,H). Then, we can 
construct level curves of R m (u,i) which are in fact the energy level curves. 
We construct level curves of R m (u,i) between two extreme values of the 
inclination: 



arccos 



H 



ft* 



"(2 



-1/2. 



< i < arccos [H] 



(13) 



being the lower one the minimum inclination to avoid the region correspond- 
ing to perihelion lower than a limit value qu m ~ a at, where the model cannot 
be applied and the higher one is due to the condition e 2 > 0. These energy 
level curves are qualitatively equivalent to the K(u, q) showed in Section [2] 
and describe very approximately the secular dynamics in the far outer solar 
system, say for a > 60 AU. At Fig. [7] we present an example of energy level 
curves calculated with this model, which corresponds with one of the maps 
of Fig. El showing the oscillations around the critical inclination in a good 
match in spite of being two quite different procedures. For a < 60 AU the 
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model reproduces qualitatively well the energy level diagrams showed in Fig. 
[5] with small differences in i c and H. 

The values of dQ/dt and dw/dt obtained from equations above are very 
small compared with the fundamental frequencies of the solar system, then 
no secular resonances should appear for any eccentr icity or any incl i natio n 
in the far TNR, which is an already know result since iKnezevic et al.l (Il99ll ) . 
But there is a particular situation of dw/dt ~ which occurs for i ~ 46° 
and % ~ 107°, where between these values the longitude of the perihelion 
has retrograde motion and prograde otherwise. This is in good agreement 
with iGallardol (j2006al ) where it was found that dw/dt ~ at i ~ 45° for the 
particular case of TNOs with a ~ 149 AU. 

More precisely, i c ~ 63° and i ~ 46° are values obtained from Eqs. (11 II) 
and ffl~2l) ignoring the higher order terms depending on E. If we take into 
account that terms these inclinations will show a dependence mainly with 
a, as showed in Fig. [HI where we also plot the TNOs with q > 30 AU and 
i > 30° which are near these inclinations. The object 2005 NU125 is near to 
the critical inclination but it is not captured around the equilibrium point, 
it has doj/dt > 0. The objects 2004 XRi 90 and 2004 DF 77 are close to the 
curve that annulate dw/dt but not enough to have their perihelion oscillating 
around a fixed value, on the contrary they have dw/dt < 0. The object 2007 
TC434 has dw/dt > 0. Finally, the objects 2006 QRiso and Eris are almost 
exactly at the curve that annulate dw/dt according to our analytical model. 
The first one exhibits large amplitude chaotic oscillations of its w, but in the 
case of Eris, its w evolve oscillating quite regularly as is showed in Fig. [9j 
From the numerical integration of Eris we obtain a circulation period of its Q 
and lj of 25.8 Myr while our analytical model predicts 26.0 Myr. We note that 
Eris, the most massive known TNO, is the only one well captured in this kind 
of motion, which can be considered as a resonance between the circulation 
frequency of its Q and u. An analogue resonance between t he cir culation 
frequency of its Q and u was found by iLibert and Delsatd (120121 ) for an 
inner particle perturbed by an eccentric giant planet with mutual inclination 
of ~ 35°, but in this case the resonant angle is u — Q, not w = u> + Q as in the 
case of Eris. However, the dynamics of Eris seems to be more complicated 
because it is very near to the three body resonance Xu + Xn — 10 A ~ 
involving Uranus and Neptune as showed in top panel of Fig. [9j Note that 
this type of secular evolution should take place also in extrasolar planetary 
systems with an external inclined low mass planet or even in a binary system 
with an external inclined planetary companion. 
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4. Kozai dynamics inside mean motion resonances with Neptune 

When imposing a resonant condition with Neptune the energy level di- 
agrams change completely. We can understand the transformation the dy- 
namics experiences when a MMR with Ne ptune is installed using a classical 



disturbing function as in iGallardol (12006a ). Considering the secular and res- 



onant terms due to Neptune in the disturbing function of a particle in the 
TNR, the equation for the time evolution of the eccentricity, for example, 
takes the form 

de 

— = / 2 sin(2w) + . . . + /o sin(er) + f[ sin(cr) cos(2w) — f[ cos(cr) sin(2w) + . . . 

1 (14) 
where a is the principal critical angle of the resonance and the / and /' are 
coefficients depending on a,e,i. There is a similar expression for di/dt. The 
critical angle is defined as 

° — (P + q)-^ - pA - qzu (15) 

where Aat,A are the mean longitudes of Neptune and the particle and the 
integers q, p are the order and degree of the resonance respectively. In the 
case of non resonant motion, the angle a is a fast circulating angle, then 
in mean it does not contribute to the time evolution of e and i making all 
terms depending on the /' to vanish and the only remaining contribution is 
due to / 2 sin(2u;), the same we obtained with our secular model of Section 
El where a; is a slow circulating or oscillating angle. But, in case of resonant 
motion, the angle a should be considered as a fast oscillating angle, and its 
contribution to the time evolution of e and i depends on the libration center, 
(To, and its amplitude. This makes an important point, as we explain below. 

From the Eq. (|T4"|) for de/dt and its analogue for di/dt it is possible to 
deduce that there should be a different behavior between a pure secular evo- 
lution (terms with /' disappear) and secular plus resonant evolution (terms 
with /' show up). Moreover, there should be a difference between exterior 
MMRs of the type 1:N and other MMRs because in the resonances of the 
type 1:N the libration center of the critical angle is cxq ^ 0°, 180°, and that 
means the terms depending on sin(cx) in mean do not vanish whereas in the 
case of the other resonances they do. Summarizing, the principal terms for 
each type of motion are the following. For secular evolution: 

^ = / 2 sim» (16) 
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for MMRs with a = 0°, 180° 



de 

-n = {h-fi< cos(a) >) sin(2w) 
at 



(17) 



and for asymmetric librations (o"o 7^ 0°, 180°, resonances 1:N): 
de 

— = (f 2 - f'i < cos(cr) >) sin(2w) + f Q < sin(cr) > +/{ < sin(cr) > cos(2o;) 

(XL 

(18) 

where < cos(a) > and < sin(er) > refer to the mean values over a libration 
cycle, which vanish for large amplitude librations converting both resonant 
cases in the secular one. The two first equations give equilibrium points at 
lo = k 90° as is usual in secular motion. But from Eq. (JIB"]) we conclude that 
the equilibrium points for resonances 1:N should be shifted from u = k90°. 

Energy level curves for a particle under the perturbation of the giant plan- 
ets and assuming a resonant motion with Neptune can be obtained by the 



Gomes et al. ( 


2005 


) and 


Gomes ( 


2011) 



numerical calculation of the energy, K, we impose a link between A at and A 
trough Eq. ( IT5l) according to an assumed time evolution of the resonant crit- 
ical a ngle, a(t), which c an be deduced from the resonant perturbing function 
as in iGallardol (j2006bl see also http://www.fisica.edu.uy/~gallardo/atlas) 
or by numerical integrations. We can calculate new energy level curves im- 
posing the resonant condition obtaining the evolution in the space (00, i) or 
(u, q) and, as we have explained above, they will be very different from 
the secular non-resonant ones analyzed in Section [21 allowing very differ- 
ent variation r anges for q and i. There are some examples in the literature 
(j Gomes et all 120051 . 120081 ; iGomesl . 1201 ll ) of energy level curves of fictitious 
particles captured in symmetric resonances showing the equilibrium points 
at u> = k90°. We also have obtained several ones for different resonances 
and libration amplitudes that confirm that behavior. To check the predicted 
behavior of the asymmetric resonances 1:N we constructed several energy 
level curves for different MMRs. For example, Fig. [TOl shows the evolution in 
the plane (u, q) of a high inclination particle captured in 1:2 resonance with 
Neptune, librating with amplitude ~ 20° around a libration center located at 
do ~ 295°, obtained from a numerical integration including all the planets. 
In the background are the energy level curves corresponding to a particle 
with the same H = 0.613 and imposing the same resonant condition that 
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the particle shows. The numerical resonant model reproduces very well the 
exact numerical integration, with small departures mainly due to variations 
in the libration amplitude that the model does not take into account. Note 
how the symmetry respect to to — k90° that the secular energy level maps 
exhibit is clearly broken. 

Figure [TT] corresponds to almost identical conditions as in Fig. [2] but im- 
posing a resonant condition corresponding to resonance 5:11 (a-d) located at 
a = 50.95 AU and 1:6 (e-h) at a = 99.45 AU. Several features can be noted 
in this figure. First, the resonant energy level maps change drastically with 
respect to the non resonant case, especially allowing large perihelion varia- 
tions not showed in the non resonant case. Second, the symmetry respect to 
to = k90° is broken in the case of resonances 1:N. This asymmetric be havior 
of the energy level curves is similar to the one obtained by first time by iKozai 
( 119851 ) in his study of Jupiter's trojans. Third, there are connections between 
perihelion outside Neptune with perihelion inside Neptune's orbit, breaking 
Neptune's barrier of the non resonant case. Fourth, the islands corresponding 
to the critical inclination are much wider and shifted from 63°. 



It has been shown by means of numerical integrations dGladman et al. 



20021 : iGomes et all 120051 iGallardol . l2006al: iGomes et all l2008h an d bv the 



const ruction of some energy level curves (IGomes et all 120051 . 120081 ; iGomesl . 



201 ll ) that the MMR+KR allows high amplitude variations of the perihelion 



distance of SDOs only for inclinations above Pluto-like inclinations. In the 
present work we can conclude that almost all substantial perihelion variations 
we observed in our energy level diagrams corresponding to a MMR+KR 
dynamics occur when the orbital inclinations are above the minimum i m ~ 
15°. Moreover, this minimum inclination necessary to connect low perihelion 
regions with high perihelion regions is greater for larger a and it occurs for 
lower values of H. For example, in Fig. UTTf h) it is shown that for resonance 
1:6 located at a = 99.45 AU, i m is approximately 15° with H = 0.7 and 
in Fig. [12] the values that can be deduced from the plots are i m ~ 20° for 
resonance 1:17 at a = 199.1 AU with H = 0.5 and i m ~ 24° for resonance 
1:30 at a = 290.8 AU with H = 0.4. 

Summarizing, when imposing the resonance condition the energy level 
curves can cross the region q ~ a at connecting by constant energy level 
curves the region q > a>N and q < allowing the biggest perihelion vari- 
ations without suffering a close encounter with Neptune (Figs. ITUl and [TT]) . 
In particular, the MMR+KR dynamics can explain the origin of some ob- 
jects with q < <2tv starting from q > without invoking close encounters 



14 



with Neptune. A pure secular, non resonant, evolution cannot generate such 
orbital transformations. Resonances 1:N which are the stronger ones in the 
TNR exhibit asymmetric energy level maps. 



5. Kozai dynamics under diffusion in semimajor axis 

It has been shown that TNOs with perihelion distances q < 36 AU experi- 
ence a diffusive process or a frankly random walk in semimajor axis due to the 
perturbations by the giant planets or by the mechanis m of resonance sticking 



contr i buting to the Scattered Disk Object population (IDuncan and Levison . 



1997 ; Gladman et al. , 20021 ) or even to the Oort cloud ( Fernandez et al. 



20041 ) . T hat upper limit for q is pushed to larger values for objects with 



larger a flGladman et all . l2002h . We re-analyzed this diffusive process in 



a wide range of q, a and % by means of numerical integrations of fictitious 
particles under the gravitational effect of the giant planets. We took 1000 
particles uniformly distributed in the phase space defined by 0° < i < 70°, 
0° < D, < 360°, 0° < u < 360°, 50 < a < 500 AU and 31 < q < f(a) where 
f(a) = 40 AU if a < 150 AU and f{a) = 30.0 + 0.085(a - 30.0 ) AU if a > 150 



AU w hich is the region susceptible of diffusion according to I Gladman et al 



(120021 ) and we followed them by 1 Gyr. 

We found the diffusion process in a is mostly restricted to the region 
q < a/27.3 + 33.3 AU and in a lesser extent it is also depending on i: small 
inclination orbits diffuse in a somehow wider region of q than large inclination 
orbits. But we found the orbital inclination is determinant in defining the 
time evolution of q: almost all particles with initial i < 20° in the region we 
studied do not experience relevant perihelion variations. In Fig. [13] we show 
in the plane (a, q) the superposition of the orbital states of all particles with 
initial i < 20°. Except for the very particular case with a ~ 100 AU there 
are no perihelion variations. On the contrary, the population with initial 
i > 20° shows evident variations in q and this is illustrated in the Fig. [TJ 
which was generated in same way as Fig. [13j For low i, paths are almost 
horizontal with no variations in q. Meanwhile, for higher i the diffusion in 
a is in general alternated with a migration in q keeping a constant, which 
indicates a capture in MMR+KR. 

We analyzed the captures in all resonances of the type 1:N up to 1:30 
and 2:N up to 2:61 (a ~ 294 AU) and we conclude that the capture in 
MMR is almost independent of the orbital inclination. In our numerical 
integrations we detected approximately 210 captures in MMR 1:N and 2:N, 
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with approximately 62% of the captures corresponding to 1:N resonances 
and 38% corresponding to 2:N resonances which is a ratio in agree ment with 
the fa c t that resonances 1:N are stronger than resonances 2:N (IGallardo . 



2006bl ; iLykawka and Mukail . 120071 ). Once a capture in MMR occurs it stops 



the diffusion in semimajor axis and a secular Kozai dynamics is installed 
inside the MMR. That secular dynamics does not allow relevant variations 
in the perihelion distances except if the inclination is high enough as we have 
explained in Section HI 

In Fig. [15] we show initial (a,i) of the particles that experienced Aq > 5 
AU due to MMR+KR, being discarded the cases of close encounters with the 
planets and also the particles ejected from the solar system. For larger a the 
required inclination necessary to generate relevant perihelion variations by 
means of MMR+KR is larger. The values of the minimum inclination that 
we deduced from the energy level maps in Section H] are in agreement with 
the numerical results given by Fig. [15] 

In the region under diffusion regime in principle neither a nor H are con- 
served and consequently the methods of Sections [2] and H] nor the analytical 
approach of Section [3] can be applied. Analyzing the variations AH in our 
numerical integrations we found that its conservation depends mostly on the 
variations Aa suffered by the particles (Fig. [16]). But, the diffusion process 
can be stopped by the temporary capture in a MMR, in particular the ones 
with enough stre ng th and width, as resonances of the type 1:N or 2:N are 



flGallardol . l2006al jbl: iLvkawka and Mukail . 120071 ; iGomes et all . 120081 ). Under 



the MMR regime the semimajor axis remains constant and the models can 
be applied. In other words, in the diffusion region the models can be applied 
during the time intervals where the motion is dominated by MMRs. 

Finally, we remark that our numerical integrations confirm that dco/dt ~ 
at i c ~ 63° and dw jdt ~ at % ~ 46°, as predicted by the analytical model. 



6. Conclusions 

We provide a series of energy level maps and analytical results describing 
the possible evolutive orbital paths for resonant and non resonant TNOs 
driven by the Kozai dynamics from Neptune and beyond. In the case of non 
resonant orbits beyond Neptune the KD does not generate relevant orbital 
variations except for orbits with inclination near the critical value i c ~ 63° 
where variations of the order of 10 AU in q are allowed and u oscillates around 
90° or 270° for a > 52 AU and around 0° or 180° for a < 52 AU. These 
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secular orbital variations always preserve the perihelion outside Neptune's 
orbit. We showed this secular KR is analogue to the mechanism that makes 
uj to oscillate in a satellite orbiting an oblate planet. Our analytical secular 
model also predicts that TNOs with i < 46° or i > 107° have a prograde 
motion of their w while orbits with 46° < % < 107° have retrograde evolution 
of their w. TNOs with orbital inclination near that limiting value should 
oscillate their longitude of the perihelion, which is the case of Eris. The 
critical inclination and the inclination that annulate the time evolution of w 
have a dependence with a, especially for a < 60 AU. 

In the case of resonant orbits, MMR+KR can generate large perihelion 
variations if their orbital inclination is greater than a minimum value i m 
that depends on a: larger a require larger i m . Contrary to the non resonant 
case, these orbital variations allow the perihelion to cross Neptune's orbit 
connecting the exterior planetary region to the interior by a secular regular 
non-encountering dynamics. 

The equilibrium points for the regime of MMR+KR are located at u — 
k90° and the energy level curves in the plane (u, q) or (u,i) are symmetric 
with respect to these values of oo as in the non resonant case, except for 
the case of resonances of the type 1:N where the symmetry is broken and 
the equilibrium points are shifted from the canonical values. This lost of 
symmetry is due to the fact that for resonances 1:N the libration center is 
not located at a = 0°, 180° as in all others MMRs. 

Particles with q < a/27.3 + 33.3 AU are inside a diffusion region in a 
where eventually could be halted by a capture in MMR. If the diffusion in a 
is greater than a few AUs, then H can not be held constant and the particle 
will not follow trajectories of equal energy. But, once the object is captured in 
a MMR the KD is installed and the evolutive paths follow curves of constant 
energy. 

Starting with a primordial population with low inclination orbits it is not 
possible to reach the present distribution of orbital inclinations of the TN 
population by means of the Kozai dynamics itself, being resonant or non 
resonant. At least, a previous excitement mechanism is necessary to enhance 
the primordial inclinations allowing the Kozai dynamics to generate a more 
diverse distribution of orbital inclinations and eccentricities in the TNR. 
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Figure 1: Trajectories of constant H (labeled lines) and the known population of objects 
with q > 36 AU (open circles). All these objects could increase or decrease their q and 
i by Kozai dynamics following curves of constant H. Data corresponding to 986 objects 
from JPL, by January 2012. 
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Figure 2: Energy level curves for a = 50 and 100 AU for different values of H showing the 
most relevant perihelion and orbital inclination variations near Neptune and beyond, q 
and i are correlated by H constant and the figures are symmetric with respect to co = k 90°. 
Note the equilibrium points located at q < 30 AU and to = 0°, 180°. Inclination in degrees. 
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Figure 3: Same as Fig. [2] for a = 200 and 300 AU. A new equilibrium point appears at 
to = 90° with a corresponding inclination of ~ 63°. For orbits with H > 0.5 no relevant 
orbital variations are allowed in this range of q. 
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Figure 4: A typical energy level map in the region q < 30 AU. Note the asymmetrical 
equilibrium islands at w ~ 70° and 110°, between the collision curves with Neptune, 
Uranus and Saturn. 
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H=0.49, a=36 AU 



H=0.4S, a=41 AU 




Figure 5: Metamorphosis of the Kozai resonance for i ~ 61° — 63°, the critical inclination, 
outside Neptune's orbit. At a = 36 AU all trajectories with q > 30 AU and H = 0.49 end 
collapsing in the collision curve with Neptune. From a = 41 AU there a ppear equilibrium 
islands at lo = 0°, 180° similar to the ones obtained bv iKuchner et al" ( 2002 '). At a — 52 
AU a new equilibrium point appear at uj = 90°, 270°. At a = 54 AU the equilibrium 
points at u) = 0°,180° disappear and fron2(Jhen on only survive the equilibrium point 
at oj = 90°, 270°. This Kozai resonance generates variations of ~ 10 AU in q and it is 
confined to i ~ 61° - 63°. 



H=0.367, a=100.4 AU 

j i 1 1 63 . 9 




30 1 1 1 1 1 1 1 59.3 

30 60 90 120 150 180 

w 

Figure 6: Example of secular long term variation in perihelion distance due to Kozai 
resonance around the critical inclination (i ~ 63°). Points correspond to a numerical 
integration of a fictitious particle by 1 Gyr. 
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30 60 90 120 150 180 

Figure 7: Energy level curves corresponding to a particle with H = 0.3 and a = 300 AU, 
the same as in Fig. [3]Jf ) but in the space (ui,i) and calculated using the analytical model 
of Section [31 
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Figure 8: The critical inclination (dashed line) and the inclination which generates 
dw/dt = (solid line) as a function of the semimajor axis, according to the analytical 
model. The TNOs with q > 30 AU are also plotted. 
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Figure 9: Orbital evolution of (136199) Eris. The orbital inclination of Eris is inside a 
range of orbital inclinations where the longitude of the perihelion oscillates as predicted 
by our analytical model. Eris is also inside the three body resonance \jj + Ajy — 10 A ~ 0. 
Orbital elements are referred to the invariable plane of the Solar System. 
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Figure 10: Points indicate the evolution in (ui,q) of a fictitious particle captured in res- 
onance 1:2 with Neptune with libration center at do ~ 295° and libration amplitude 20° 
obtained by numerical integration in an interval of 10 Myrs. In the background are the 
energy level curves for a particle with the same H = 0.613 and same oq and amplitude as 
the particle. 
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Figure 11: Energy level curves for orbital configurations similar to Fig. [5] but imposing 
a resonant condition. From (a) to (d) correspond to resonance 5:11 assuming <7o = 180° 
and amplitude 20°. From (e) to (h) correspond to resonance 1:6 assuming cto = 60° and 
amplitude 20°. Compare with Fig. [2j For growing libration amplitude the energy level 
curves tend to the secular non resonant case. 
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Figure 13: Superposition of the orbital states of a numerical integration by 1 Gyr of the 
outer Solar System plus 273 fictitious particles with initial 0° < i < 20°. The dashed 
line corresponds to q = a/27.3 + 33.3 AU which approximately defines the limit for the 
diffusion region in semimajor axis. Diffusion in a can be halted by a MMR but once a 
capture in MMR occurs the Kozai dynamics does not generate relevant changes in q due 
to the very low orbital inclination (with an extraordinary exception of a particle captured 
in resonance 1:6 at a <~ 100 AU). 
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Figure 14: Superposition of the orbital states of a numerical integration by 1 Gyr of 
the outer Solar System plus 727 fictitious particles with initial 20° < i < 70°. The 
diffusion process is frequently stopped by a capture in MMR+KR allowing large perihelion 
variations. However, larger a require larger i to generate relevant perihelion variations (see 
Fig. [T5|) . The dashed line is the same limiting line as in Fig. [T3] 
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Figure 15: Initial (a, i) for particles experiencing perihelion variations Aq > 5 AU due 
to Kozai resonance. Orbits with larger a require larger i to generate relevant perihelion 
variations by means of MMR+KR. Particles having close encounters with the planets were 
not considered in this plot, the inclination was taken with respect to the invariable plane 
of the Solar System and the semimajor axis is baricentric. 
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Figure 16: H fina i - H initia i versus \af ina i - a init iai\ for particles surviving all the in- 
tegration. Particles with diffusion in a up to few UAs conserve their H but otherwise 
considerable variations in H occur and the evolutive paths of these particles will not be 
confined to the trajectories of H = constant. For the calculation of H the inclination was 
considered with respect to the invariable plane of the Solar System. 
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